On Hamiltonian intermittency in equal mass three-body problem 
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We demonstrate that both kinds of the Hamiltonian intermittency exert an influence on 

o 

' the disruption statistics in the equal mass three-body problem. Studying initially-resting 

^ , triple systems we found a narrow region in the vicinity of the strong chaos, where the 

influence of the second kind Hamiltonian intermittency (T^ '^^^) trajectories cause the integral 
psj ■ distribution to distort enough to be detected. We fitted the integral distribution with both 

power-laws (T^ "^^^ and ^^^) taken into account, and found an excellent agreement between 
the fit and observed integral distribution. 
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PACS numbers: 05.45.-a, 05.45. Pq, 45.50.Pk, 95.10.Ce, 95.10.Fh 
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exponential tail in the disruption statistics of the 

equal mass three-body problem. 
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The three-body problem, although being the 
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tH- ' not have an analytical solution in the gen- , ^ ^ i j_ • i • i • i i , 

' demonstrated that m the hierarchical three- body 



^ ■ eral case, which makes it possible to gain an , , r ■ 1 t t -t^ 

. problem the decay of the survival probability 



understanding of its dynamics mostly through 
the means of numerical studies. Those studies 



is heavy-tailed with a power index equal to 
-2/3, which corresponds to the first kind of 



, , , started even before any computers were devel- tt -i^ ■ • ^ -^^ 1 j_ -i 

_ Hamiltonian intermittency; the power-law tails 

H ' oped and trace back to the beginning of the 20th . , , , .111 11 ■ . 1 ■ 

Cd ■ n 00 equal mass three- body problem with m- 

1|. With the development of computers ,. t ^ ^i ^ ^ ^^ ^ i 

' dices close to that value were recently reported 



century 

the effectiveness of these studies increased, but, 
due to the insufficient amount of data, this does 



by Orlov et al. 



(they also appear in the 



decay of the survival probability of the exited 

not led to an immediate breakthrough in the ^ Jx i 1 ■ i r j_i tt -t^ ■ 

atoms pj). ihe second kind of the Hamiltonian 

understanding of the disruption process. Val- . ^ .^^ i- ^ ■ ^ r 

jr intermittency predicts the existence of a power- 

assumed that the tail of the disruption , , -i , . , 1,0/0 ■ 1 

law tails with index equal to -o/2; a power index 

distribution should be exponential, though this 



close to this value was reported by Shevchenko 



assumption conflicts wit 
result by Agekian et al. 



1 the earlier theoretical 
that the average life- 
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and Scholl 



8| in the 3/1 Jovian resonance (Sun- 



Jupiter-Asteroid problem). Other than that. 



time of an isolated triple system should be infi- 
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it appears that power-law tails are common for 
the disruption statistics of the Hamiltonian svs- 

Q 

terns of different nature (see [9j] and references 
therein). 

As it was mentioned, both kinds of Hamilto- 
nian intermittency were seen in the three-body 
problem, hence, they both are "native" to this 
problem, yet, no one reported seing both of them 
in one set. That is why we decided to thoroughly 
investigate the equal mass three-body problem 
and demonstrate that both of them are present. 
In the next section we describe our numerical 
method and the initial data set, then we report 
the results and discuss them in the final section. 

II. NUMERICAL METHOD 

To perform our calculations we use the ver- 
sion of the Aarseth-Zare triple code [10(1. In 

n 

contrast with our previous paper we used 
a "standard" way to define the initial condi- 
tions for the initially-resting triple system, used 
in Agekian et al. [s], where the authors demon- 
strated that the approach used cover all possi- 
ble configurations of three bodies. It is demon- 
strated in Fig.[T]- that two bodies are placed in 
points A and B, while the third body could be 
placed in any point in the gray region. The gray 
region is bounded by the unit circle centered in 
A from the right, x = 0.5 from the left and x- 
axis from the bottom. Scanning over all possi- 
ble initial positions of the third body gives us a 
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FIG. 1: Initial positions of the masses: first and sec- 
ond are located in points A and B while third mass 
could be placed in any point in the gray region 

complete set of data (as we use the equal-mass 
problem, that is enough; should we use different 
masses, we have to permute the masses between 
A, B, and the "running" position, and use all 
the resulting time maps together). Initially, all 
three masses are at rest, but immediately after 
they start moving in the gravitational potential 
they govern. The calculation of their motion 
lasts until the system disrupts (at time T^) or 
T = 5000 is reached; the value for the maximal 
T is adopted from [1 1| , where we saw that by this 
time '^^'^ is clearly reached in the equal-mass 
case, our current research proved that as well. 
The disruption of the system is determined by 
the hyperbolicity (positivity of mechanical en- 
ergy of the disrupting binary) at a distance 50 
times larger than the current semi-major axis of 
the final binary. This way of defining is by de- 
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fault used in the triple code that we are using 
and it is quite wide-spread in studying the three- 
body problem in astronomy. 

III. RESULTS 

We have modelled over 30 million initial 
positions for the third mass that are evenly- 
distributed over the gray area in Fig.[Tl The re- 
sulting Fa{Td) integral distribution (the number 
of trajectories with the disruption time that is 
greater than T^) was fitted by the power-law re- 
lation Fa{Td) = AxT^ with steps AT = 10 and 
AT = 100; the result is presented in the upper 
panel of Fig.O AT = 10 as the red (dark gray 
if grayscaled) and AT = 100 as the thick black 
line. The strong noise of the first curve is caused 
by an insufficient number of data points in each 
bin (AT = 10), one can easily see that the sec- 
ond curve is much smoother and still, it has vari- 
ations. They both approach /3 = —0.7 ... — 0.75, 
which is pretty close to /3 = —2/3 that is ex- 
pected from the first type of the Hamiltonian 
intermittency. On the other hand, the mini- 
mum at low Td points to /3 = —3/2 (although, it 
does not reach it, ending at (3 ^ —1.06) that is 
expected from the second type of the Hamilto- 
nian intermittency, hence we investigate it fur- 
ther on. One note needs to be taken - with 
AT > 200 . . . 400 minimum at low T^ could be 
easily missed and we believe that large AT was 
a reason why /3 = —3/2 was not detected in this 



problem earlier. 

Then we fit the Ta(Trf) curve as a sum of both 
contributions Fa{Td) = x T^' + A2 x T^' 
and this fit gives us an excellent result: 
Ai = 3.93812 X lO'^ ± 5.297 x 10^ (1.345%), 
A2 = 1.2327 X 10^ ± 3.707 x lO'^ (3.007%), 
/3i = -0.666399 ± 0.001526 (0.229%), and 
/32 = -1.501 ± 0.008454 (0.5633%). One can see 
that the powers are almost exactly equal to the 
predicted values; in the bottom panel of Fig.[2] 
we presented the Fa{Td) curve (thick gray), 
T^ contribution (dotted green), T^ (red 
dashed), and the sum of two contributions (thin 
blue). One can see that the resulting curve fits 
Fa{Td) quite well, while T^ alone coincide 
with Fa{Td) only starting from thousands of Td- 

IV. CONCLUSIONS 

We have considered the three-body prob- 
lem with unit masses that are initially at rest 
and demonstrated that both terms - T^ "^^^ and 
T^ ^'^^ - contribute into the statistics of the dis- 

—3/2 

ruption. The reason why T^ contribution re- 
mained undetected in the equal-mass problem is 
because the corresponding trajectories are stuck 
to the strong chaos region at low Td- The chaotic 
pattern ends at the low values of T^, so the tra- 
jectories that are stuck to its border manifest 
themselves strongest at low Td as well (though, 
not solely, since it is owing to the influence of 

—3/2 —2/3 

T^ contribution that the power of T^ is 
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FIG. 2: Upper panel: the power index of the Fa{Td) power-law fit {Fa{Td) (x T^)- The red (dark gray if 
grayscaled) curve corresponds to AT = 10 step while the thick black curve - to AT — 100. Bottom panel: 
best- fit of the Fa{Td) curve with both contributions: the thick gray line - Fa{Td), the green dotted line - 
contribution, the red dashed line - , the thin blue line - the sum of two contributions. See the 



text for exact values of powers. 
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not found exactly). Fitting the FaiTd) curve 
clearly demonstrated the presence of both the 
contributions (T^ and ) , and the powers 
were confirmed with great precision. We believe 
that the reason why the power indices found in 
the earlier papers on this subject (e.g.,jiQi) were 
somehow lower than the expected —2/3 was the 
interference from the "^^"^ term. Also, high val- 
ues for the maximal integration time usually lead 
to high values of AT - the size of bin to plot 
the distribution. If AT > 200 . . . 400, then the 
whole deep in upper panel of Fig. [2] would be in- 



side one bin, and so it could be undetected - we 
believe that this was, at least, partially, a rea- 
son why both indices simultaneously remained 
undetected. 
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